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We present rigorous performance bounds for the optimal dynamical decoupling pulse sequence 
protecting a quantum bit (qubit) against pure dephasing. Our bounds apply under the assumption of 
instantaneous pulses and of bounded perturbing environment and qubit- environment Hamiltonians. 
We show that if the total sequence time is fixed the optimal sequence can be used to make the 
distance between the protected and unperturbed qubit states arbitrarily small in the number of 
applied pulses. If, on the other hand, the minimum pulse interval is fixed and the total sequence 
time is allowed to scale with the number of pulses, then longer sequences need not always be 
advantageous. The rigorous bound may serve as testbed for approximate treatments of optimal 
decoupling in bounded models of decoherence. 

PACS numbers: 03.67.Pp, 82.56.Jn, 76.60.Lz, 03.65.Yz 



I. INTRODUCTION 

Quantum systems tend to rapidly decohere due to the 
coupling to their environments, a process which is es- 
pecially detrimental to quantum information processing 
and high resolution spectroscopy [l| . Of the many meth- 
ods which have been proposed in recent years to over- 
come the damage caused by decoherence, we focus here 
on dynamical decoupling (DD), a method for suppress- 
ing decoherence whose origins can be traced to the Hahn 
spin echo In DD, one applies a series of strong and 
frequent pulses to a system, designed to decouple it from 
its environment [3|-|6|. Recently, it was discovered how 
to optimally suppress decoherence of a single qubit us- 
ingDD, under the idealization of instantaneous pulses 
041] ■ One of us found an optimal pulse sequence (later 
dubbed Uhrig DD, or UDD) for suppressing pure dephas- 
ing (single-axis decoherence) of a qubit coupled to a bo- 
son bath with a hard frequency cut-off Q. In UDD one 
applies a series of N instantaneous 7r pulses at instants 
tj (j € {1, 2 . . . N}), with the instants given by tj = Tdj 
where T is the total time of the sequence and 



5j =sin 2 (j7r/(27V + 2)). 



(1) 



By optimal it is meant that with each additional pulse 
the sequence suppresses dephasing in one additional or- 
der in an expansion in T, i.e., N pulses reduce dephasing 
to 0(T N+1 ). The existence and convergence of an expan- 
sion in powers of T, at least as an asymptotic expansion, 
is a necessary assumption [Sj, [10| . 

The UDD seq uence was first conjectured [ill, EH and 
then proven [l(| to be universal, in the sense that it ap- 
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plies to any bath causing pure dephasing of a qubit, not 
just bosonic baths. The performance of the UDD se- 
quence was tested, and its advantages over other pulse 
sequences confirmed under appropriate circumstances, in 
a series of recent experiments |13l - tl5| . Its limitations as 
a function of sharpness of the bath spectral density high 
frequency cut-off (l6l | and as a function of timing con- 
straints [13] have also been discussed. 

In order to suppress general (three-axis) decoherence 
on a qubit to all orders concatenated sequences are 
needed [HI, [l!| , whose efficiency can be improved by us- 
ing UDD building blocks [20j | . A near optimum sup- 
pression is achieved by quadratic dynamic decoupling 
(dubbed QDD). This scheme was proposed and numer- 
ically tested in Ref. [8( and analytically corroborated in 
Ref. [9]. In QDD. a sequence of (JV + l) 2 pulses compris- 
ing two nested UDD sequences suppresses general qubit 
decoherence to 0(T N+1 ), which is known from brute- 
force symbolic algebra solutions for small N to be near- 
optimal Q. 

While rigorous performance bounds have been derived 
previousl y for p eriodic and concatenated DD pulse se- 
quences [2ll423||. no such performance bounds have yet 
been derived for optimal decoupling pulse sequences, in 
particular UDD and QDD. In this work we focus on 
UDD and obtain rigorous performance bounds. We post- 
pone the problem of finding rigorous QDD performance 
bounds to a future publication. Our main result here 
is an analytical upper bound for the distance between 
UDD-protected states subject to pure dephasing and un- 
perturbed states, as a function of the natural dimension- 
less parameters of the problem, namely the total evo- 
lution time T measured in units of the maximal intra- 
bath energy Jo, and in units of the system-bath coupling 
strength J z . The bound shows that this distance (tech- 
nically, the trace- norm distance), can be made arbitrar- 
ily small as a function of the number of pulses N, as 
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(1/M)(J T + J^T)^. This presumes that the bounds 
J a (a € {0, z}) are finite, an assumption which will fail 
for unbounded baths, such as oscillator baths. In such 
cases, which includes the ubiquituous spin-boson model, 
our analysis does not apply. Alternative approaches, such 
as those based on correlation function bounds [23(, are 
then required. 

We begin by introducing the model of pure dephasing 
in the presence of instantaneous DD pulses in Section [TT1 
We then derive a general time-evolution bound in Section 
IIII1 without reference to any particular pulse sequence. 
In Section IIVI we specialize this bound to the UDD se- 
quence. Then, in Section |Vj we obtain our main result: 
an upper bound on the trace-norm distance between the 
UDD-protected and unperturbed states. In Section IVII 
we analyze the implications of this bound in the more 
realistic setting when only a certain minimal interval be- 
tween consecutive pulses can be attained. Certain techni- 
cal details are presented in the Appendix, including the 
first complete universality proof of the UDD sequence, 
which does not rely on the interaction picture. 



II. MODEL 

We start from the general, uncontrolled, time- 
independent system-bath Hamiltonian for pure dephas- 



ing 



ffimc = Is ® B + cr z (g) B z , 



(2) 



where Bo and B z are bounded but otherwise arbitrary 
operators acting on the bath Hilbert space Hb , Is is the 
identity operator on the system Hilbert space Us, and 
(7 z is the diagonal Pauli matrix. The bath operator B z 
need not be traceless, i.e., we allow for the possibility of a 
pure-system term a z CS> Ib in the system-bath interaction 
term H^b '■= cr z <£> B z . Such internal system dynamics 
will be removed by the DD pulse sequence we shall add 
next, along with the coupling to the bath. However, the 
assumption of pure dephasing means that we assume that 
the level splitting of the system, i.e., any term propor- 
tional to a± (with cr± being cos(ip)a x + sin(</?) cr y for arbi- 
trary tp) acting on the system, is fully controllable. Oth- 
erwise the model is one of general decoherence, and our 
methods require a modification along the lines of Refs. [8[ 
and @. If the system described by Eq. @ is subject to 
N instantaneous ir pulses at the instants {tj := T8j}^ =1 
about a spin axis perpendicular to the z-axis, i.e., if 
the Hamiltonian H^B-it) = f YljLi — %) (T -L ® is 
added to H unc , the interaction picture ( "toggling- frame" ) 
Hamiltonian H tog (t) = Ul, D (t)H unc UBB(t) reads 



H tog (t) = I s ® B Q + f(t)a z ® B z , 



(3) 



where the unitary J7dd(0 alternates between 7s (8>/b and 
ct_l ® Ib at the instants {tj}^L 1: and consequently the 
"switching function" fit) = ±1 changes sign at the same 
instants. 



We shall also need the magnitudes of the two parts of 
the Hamiltonian 



Jo := ||.Bn|| < oo, J z := \\B Z \\ < oo, 



(4) 



where |j • |j is the sup-operator norm (see Appendix |A| . 
There are certainly situations where either J or J z can 
be divergent (e.g., Jo in the case of oscillator baths). In 
such cases our bounds will not apply, and other methods 
such as correlation function bounds are more appropriate 
(see, e.g., Ref. [23| for such bounds applied to DD). 



III. TIME EVOLUTION BOUNDS 

We aim to bound certain parts of the time evolution 
operator induced by Htog(t) 



U{T)=Tenp(^-i H tos (t)dt 



(5) 



where T is the time-ordering operator. 

Standard time dependent perturbation theory provides 
the following Dyson series for U(T) 



U(T) 



n=0 



-iT) n 



E 

{a;dim(a)— n\ 



(6a) 
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'o Jo 
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ds 2 f a2 (s 2 ) \ dsif ai (si) (6b) 
10 Jo 
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(6c) 



where dim(a) = n is the dimension of the vector a. 
The identity Is in the Hilbert space of the qubit/spin 
is denoted by cto- In all sums over the vectors a their 
components ctj take the values or z. In this way, 
the summation includes all possible sequences of Bq and 
B z . The function /o(s) is constant and equal to 1 while 
f z (s) := f(sT) takes the valus ±1. We use the dimen- 
sionless relative time s := t/T so that all dependence on 
T appears as a power in the prefactor. Note that the 
coefficients F3 do not depend on T. 

In order to find an upper bound on each term F3 Qs 
separately we proceed in two steps. First, we use \f a \ = 1 
to obtain 



\F 3 \ < / ds n / ds„_i 
Jo Jo 
1 



S3 PS2 

ds2 / ds\ 

n Jo 



Second, we use Eq. ^ and \\cr a \\ — 1 to arrive at 



\Qs\\ < l[ja 



n-k(u) jk(S) 



(7) 



(8) 
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where we used the submultiplicativity of the sup-operator 
norm (see Appendix lAl . The number k(a) stands for 
the number of times that the factor J z occurs. Standard 
combinatorics of binomial coefficients tells us that the 
term J^~ k J z occurs n\/{k\{n— k)\) times in the sum over 
all the vectors a of given dimensionality n in (|6a|) . Hence 
each term of the time expansion of U(T) is bounded by 



{a;dim(a)— n} 



1 



k\(n 

k=0 y 



7" 

k)\ J ° 



-k T k 



(9) 



We therefore define the bounding series 



S(Jo, Jz 



1 



n=0 k=0 

exp((J + J Z )T) 



k\{n-k)\ 



jn — k jk 
J J z 



(10a) 
(10b) 



It then follows from Eq. ^ that each multinomial in Jo 
and J z of the expansion of S(Jq, J z ) is an upper bound 
on the norm of the sum of the corresponding multinomial 
in the operators Bq and B z of the expansion of U(T) in 
Eq. (|6a[) . This is the property which we will use in the 
sequel. 



IV. BOUNDS FOR DEPHASING 

From <t z — Ig it is obvious that only the odd powers in 
B z contribute to dephasing while the even ones do not. 
Hence we split U(T) as 



U(T) = Is <8 B+(T) + g z ® B-(T) 



(11) 



where the operators B± act only on the bath while Ig and 
a z act only on the qubit. The operator B+ comprises all 
the terms with an even number k of a z ® B z , i.e., with an 
even number of J z in the bounding series S(Jq, J z ). The 
operator _B_ comprises all the terms with odd number k 
ola z ®B Zl i.e., with an odd number of J z in the bounding 
series S(Jo, J z ). Hence to bound the time series of B_(T) 
term by term we need the the time series of the odd part 
of S(J , J z ) in J z . This, from (flOE) is: 



S-(J , Jz) = exp(J T) sinh( J Z T). 



(12) 



The time series of S~(Jq, J z ) provides a bounding series 
of B- (T) term by term. Hence we define 



1 d k 

dk := ~k\dT k ~ S ^ J °' Jz ^ 



T=0 



(13) 



such that S-(J , Jz) = EZo d kT k - 

We know from the proof of Yang and Liu [10( that in 
the _Bo-interaction picture a UDD sequence with TV pulses 
(which we denote by UDD (TV)) should make the first TV 
powers in T of BJT) vanish, i.e., B_(T) = 0{T N+1 ). 
However, since the Yang-Liu proof does not directly ap- 
ply to our discussion, we provide a complete version of 



this proof which avoids the i?o-interaction picture in Ap- 
pendix [Bj The remaining powers are bounded by the 
corresponding coefficients dk of S- . Thus the expression 



A N := J2 d * Tk 



(14) 



k=N+l 



provides an upper bound for B_(T) if UDD(iV) is ap- 
plied: 



\B-{T)\\ < A 



JV- 



(15) 



Due to the obvious analyticity in the variable T of 
S-(Jq,J z ) as defined in (TT21 we know that the residual 
term vanishes for N — > oo, i.e., 



lim A 



N 



0. 



(16) 



This statement holds true irrespectively of the values of 
Jo and J z , as long as they are finite. 

We can obtain a more explicit expression for A^r. Be- 
sides the dimensionless number of pulses N the bound 
Ajv depends on JoT and on J Z T. It is convenient to 
introduce the dimensionless parameters 

e := J T, r, := J Z /J Q (17) 

instead. In terms of these parameters we have 

S- (77, e) = exp(e) sinh(e?7). (18) 

From the series 



cxp(e) sinh(ery) = i[e £< - 1+ ''' — e 



(19a) 

E5?[( 1 +^-( 1 -^ ( 19b ) 

1=0 

00 

(19c) 



00 / 



1=0 



with 



we obtain 



PI(V) :=^j [(1 +»?)'- (l-»7)'] 



(20) 



A N ( v ,e) = P*(v)e n (21a) 

n=N+l 

- PN+1 { V )e N+l +0{e N+2 ). (21b) 

This, together with the bound ([T5|). is our key result: 
it captures how the "error" ||S_(T)|| is suppressed as a 
function of the relevant dimensionless parameters of the 
problem, ij, e, and N. Note that convergence for N — > 00 
is always ensured by the factorial in the denominator, 
irrespectively of the values of e and rj as long as these are 
finite. 
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For practical purposes it is advantageous not to com- 
pute A N by the infinite series in (|21al) , but by 



A N foe) = S- (r),e) 



N 

£ 

n=0 



(22) 



which can easily be computed by computer algebra pro- 
grams. Figures [1] and [2] depict the results of this compu- 
tation. Consider first Fig.Q] Each curve shows Ajv(?7,e) 
as a function of e, at fixed r\ and N. The error ||S_(T)|| 
always lies under the corresponding curve. Clearly, the 
bound becomes tighter as e decreases. Moreover, the 
more pulses are applied (the different panels) the higher 
the power in s and thus the steeper the curve. Addition- 
ally, the curves are shifted to the right as N increases. 
Clearly, then, a larger number of pulses improves the er- 
ror bound significantly, at fixed e and 77. This effect is 
even more conspicuous in Fig. [5J where 77 is fixed in each 
of the two panels, and the different curves correspond 
to different values of N. The vertical line intersects the 
bounding function at progressively lower points as N is 



3 0(2!"] 




FIG. 1: (Color online) The bounding function Ajv as a func- 
tion of e = JqT, as given in Eq. (|14|l for various numbers of 
pulses N and various values of the parameter rj = J z /Jq € 
{0.01, 0.1, 1, 10, 100}, with r\ increasing from the rightmost 
curve to the leftmost curve in each panel. 



V. DISTANCE BOUND 

Intuitively, we expect the bound on ||_B_(T)|| derived 
in the previous section to be sufficient to bound the effect 
of dephasing. However, to make this rigorous we need a 
bound on the trace-norm distance D[p§(T), Pg(T)] be- 
tween the "actual qubit" state 



p s (T) :=ti B \p SB (T)] 
and the "ideal qubit" state 

p° s (T) :=tr B [p° SB (T)] 



(23) 



(24) 
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FIG. 2: (Color online) The bounding function Am as a func- 
tion of s = JqT, as given in Eq. (|21a|) for various numbers of 
pulses N £ {2,5, 10,20}, at fixed values of 77. In each panel 
the curves become steeper as N increases. 

where p^ B (T) is the time-evolved state without coupling 
between qubit and bath. The partial trace over the bath 
degrees of freedom is a map from the joint system-bath 
Hilbert space to the system-only Hilbert space (see Ap- 
pendix [A}, and is denoted by trs- As we shall see, the 
term 1$ <8> B + (T) in Eq. (fTTj) indeed has a small, and in 
fact essentially negligible effect. 

To obtain the desired distance bound we consider a fac- 
torized initial state Pg B (0) = \tp)(ip\ <8> Pb, which evolves 
to Psb(T) = U(T)p^ B {0)W{T) when the system-bath 
interaction is on (the "actual" state), or to /Og B (T) = 
Is ® U B {T)pl B (0)I s ® U B (T) when the interaction is off 
(the "ideal" state). The unitary time evolution operator 
without coupling reads 

U B (T) := exp(-tTBo), (25) 

where Bq is the pure-bath term in Eq. ([2]). The initial 
bath state pb is arbitrary (e.g., a mixed thermal equilib- 
rium state), while the initial system state is pure. Let us 
define the correlation functions 



(T) :=tr \B a {T)p B BUT) 



(26) 



where a, /3 G {+,—}, and where all operators under the 
trace act only on the bath Hilbert space. Explicit com- 
putation (see Appendix [Ct then yields: 

D[p s (T),p° s (T)] (27) 

< \{\H + {T) - 1| + |6 + _(T)| + \b. + (T)\ + |6__(T)|). 

We will show that 6++ is very close to 1 while the 
other b a p quantities are small in the sense that they are 
bounded by Eq. (fT5)) . 

First note from the unitarity of Eq. (fTT|) that 

I = U^U (28) 

= J s <g> (B f + B + + + a z ® {bIb+ + B\B-) 
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where we omitted the time dependence T to lighten the 
notation. Hence we have 

I = B\b+ +BiB- (29a) 
= B\B- +BlB + . (29b) 

It follows that (i\B\_B + \i) = \\B + \i)\\ 2 < 1 for all nor- 
malized states because (i\B-B^_\i) = 2 is non- 
negative. Thus in particular max|||j\|i = i ||-B_|_|i)|| < 1, and 
we can conclude that 

\\B+\\ < 1. (30) 

Cyclic invariance of the trace in b a p together with 
Eq. (|29a|) and the normalization tr[ps] = 1 immedi- 
ately yields + 6 = 1, while the combination with 

Eq. (|29b|) implies 6 + _ + = 0. Hence Eq. ([27} can be 
simplified to 

D[p s (T),p° s (T)]<\b + -(T)\ + \b-(T)\. (31) 

To obtain a bound on the correlation functions b a p we 
use the following general correlation function inequality 
(for a proof see Appendix |D|) : 

\tr[Q PB Q']\ < IIQ'IHIQII, (32) 

which holds for arbitrary bounded bath operators Q,Q'. 
Applying Eq. J32J to Eq. <(26]) yields 

|6__(T)| < ||S-(T)|| 2 , (33a) 
\b + .(T)\ < \\B + (T)\\\\B_(T)\\ (33b) 
< \\B-(T)\\, (33c) 

where in the last inequality we used Eq. ([3U|) . 

Summarizing, together with Eqs. (|15[) and (|3ip we 
have obtained the following rigorous upper bound for the 
trace-norm distance 

D[p s (T),p° s (T)) <mm[l,A N (r),e) + A 2 N (r,,e)]. (34) 

This upper bound completes our main result. Since as we 
saw in Eq. (J21b| A N (r),e) = p N+1 {r))e N+1 + 0{e N + 2 ), 
the appearance of the squared term in Eq. (j34|) (whose 

origin is \b (T)\) is not relevant in the sense that even 

in the presence of this term the bound 

D[ Ps (T),p° s (T)} < pn + i(v)s N+1 + 0(e N + 2 ) (35) 

holds. Hence the result of Eq. (|21a[) depicted in Figs. Q] 
and [JJ provides the desired result. Ignoring the A^ term 
in Eq. (fM]) . we note that Figs. Q] and [2] also reveal the 
limitations of our bound when e or r\ are too large for a 
given value of N: For any pair of states it is always the 
case that D < 1, so that as soon as An = 1 the bound 
no longer provides any useful information. 

Note further that the results shown in Fig.[T]are quali- 
tatively similar to the results obtained for the analytically 
solvable spin-boson model for pure dephasing Q . Heuris- 
tically, the necessary identification is Jo = wd where ud 



is the hard cutoff of the spectral function and rj oc a 
where a is the dimensionless coupling constant for Ohmic 
noise. We stress that the advantage of Eq. (l34|) com- 
pared to the analytically exact results in Ref. [7| is that 
it holds rigorously for a large class of pure dephasing 
models, namely those of bounded Hamiltonians. 



VI. ANALYSIS FOR FINITE MINIMUM PULSE 
INTERVAL 

So far we have essentially treated the total time T 
and the number of pulses N as independent parame- 
ters. This is possible when there is no lower limit on 
the pulse intervals. However, in reality this is never 
the case and in this section we analyze what happens 
when there is such a lower limit. Note that it follows 
from Eq. (JXJ) that the smallest pulse interval is the first: 
t\ = T sin 2 (tt/(2N + 2)). Let us assume that t\ is fixed, 
so that, given t\ and TV, the total time is 

T(N) = t iq (N), (36a) 

9 W := esc 2 (^) . (36b) 

For large N we can expand the esc 2 function to first 
order in its small argument, yielding 

q{N) =(^L±iy + 1 - + 0(N- 2 ), (37) 

which shows how the total time grows as a function of 
N at fixed minimum pulse interval t\. Along with 77, the 
relevant dimensionless parameter is now 

£1 := Joh, (38) 

instead of e = q{N)e\ . We can then rewrite the bounding 
function (|21ap in terms of these quantities as 

00 

A N (r),ei)= P n (v)q n (n)E^ (39) 

n=N+l 

Considering now the large N limit of the first term in 
this sum, we have 

p N+1 {ri)q N+1 {N + l)e^ » JL (^V) 2W [(1 + „) Bl ] w 

w {cN) N (40) 

where we kept only the leading order terms and neglected 
all additive constants relative to N, and in Eq. (l4Cfl) used 
Stirling's approximation n\ ~ (n/e) n . The constant c is 
i(^) 2 e(l + 77)^1 . We thus see clearly that for fixed t\ it 
becomes counterproductive to make N too large, since 
no matter how small c is, for large enough N the fac- 
tor TV^ will eventually dominate. This reflects the com- 
petition between the gains due to higher order pulse se- 
quences and the losses due to the increased coupling time 
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to the qubit allotted to the environment. Similar con- 
clusions, delineating regimes where increasingly long DD 
sequences become disadvantageous, have been reported 
for periodic [H and concatenated 0, [H, Hi DD 
pulse sequences, as well as for the QDD sequence Q. 

These conclusions are further illustrated in Fig. [3l 
where we plot the bound Ajv(?7, £i) by replacing s with 
e±q(N) in Eq. ()21a|) . This figure should be contrasted 
with Fig. [2j The most notable change is that increas- 
ing N now no longer uniformly improves performance. 
Whereas in Fig. [5] the curves for different values of N all 
tend to converge at high values of e, in Fig. [3] a high N 
value results in a steeper slope, but also moves the curve 
to the left. Thus, for a fixed value of ei it can be advan- 
tageous to use a small value of N (e.g., for rj = 0.01 and 
£i =0.1 the N = 2 curve provides the tightest bound). 
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case when the uncontrolled Hamiltonian ([2]) is replaced 
by H unc = Is®Hb+ct±®B, and the UDD pulse sequence 
consists of rotations about the spin-z axis. A practical 
implication is that the bounds found here can be used to 
check numerical and approximate calculations. Such cal- 
culations must obey our mathematically rigorous bounds, 
so that a testbed is provided. 

Furthermore, a number of interesting generalizations 
and extensions of our results readily suggest themselves. 
One is to consider rigorous bounds for finite pulse-width 
UDD sequences. It is already known how to construct 
such sequences with pulse-width errors which appear only 
to third order in the value of the pulse width [25J, but 
no rigorous bounds have been found. Another important 
generalization, as mentioned above, is to the QDD se- 
quence for general decoherence d, We expect that 
techniques similar to the ones we introduced here will 
apply to both of these open problems. Yet another di- 
rection, which will require different techniques, is to find 
rigorous UDD performance bounds for unbounded baths, 
such as oscillator baths. It is likely that a correlation 



function analysis similar to that performed in Ref. 23 1 
for periodic and concatenated DD sequences will prove 
useful in this case. 
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FIG. 3: (Color online) The bounding function A at as a func- 
tion of ei = ti Jo (where ti is the smallest pulse interval), at 
fixed values of rj. The number of pulses TV £ {2, 5, 10, 20} is 
varied from curve to curve. In each panel the curves become 
steeper as N increases. 



VII. CONCLUSIONS 

We have derived rigorous performance bounds for the 
UDD sequence protecting a qubit against pure dephas- 
ing. The derivation is based on the existence of finite 
bounds for the relevant parts of the Hamiltonian, cap- 
tured in the dimensionless parameters e and rj. Under 
this assumption the bounds show rigorously that dephas- 
ing is suppressed to leading order as (l/iV!)[e(l + rj)] N 
We consider it a vital step to know that irrespectively of 
any details of the bath, except for the existence of finite 
bounds, a large number N of pulses is always advanta- 
geous at fixed T - at least under the idealized assumption 
of perfect and instantaneous pulses. 

An immediate corollary of our results is that identical 
bounds apply for the case of the UDD sequence protect- 
ing a qubit against longitudinal relaxation. This is the 



Appendix A: Norms and Distances 
1. Trace and partial trace 

We deal only with linear trace-class bounded opera- 
tors that map between separable Hilbert spaces in this 
work. A Hilbert space Ti is separable if and only if it 
admits a countable orthonormal basis. A bounded lin- 
ear operator A : Ti n- Ti, where Ti is separable, is said 
to be in the trace class if for some (and hence all) or- 
thonormal bases {|fc)}fc of T-L the sum of positive terms 
5Zi.(fc|yl''' A\k) is finite. In this case, the sum 5Z fe (fc|-<4-|fc) 
is absolutely convergent and is independent of the choice 
of the orthonormal basis. This value is called the trace of 
A, denoted by tr(A). Whenever we use the symbol tr in 
this work, we mean the trace over the full Hilbert space 
the operator the trace is taken over is acting on. 

Now consider two separable Hilbert spaces Ti\ and 
Ti.2 and let A : Ti H> Ti denote a linear trace-class 
bounded operator acting on the tensor product Hilbert 
space Ti := Ti\ ®Ti%. Let {\ki)}k denote an orthonormal 
basis for Tii, where i £ 1,2. The partial trace opera- 
tion over the first (second) Hilbert space is a map from 
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H to the second (first) Hilbert space, and has the op- 
erational definition trj(A) := Ylk- (fcj|j4|A;j}. When A is 
decomposed in terms of the two orthonormal bases as 
A = Y. kuKM ^Mmmih){k' 1 V 2 \, where \k x l 2 ) := 
\k\) ® etc., the partial trace over % 2 can written as 
tr 2 (^) = £ fcl , fc; , i2 (M 2 |^|^2}|fci}(fci|. This makes it 
clear that tr 2 (A) is an operator that acts on Hi. 



2. Sup-operator norm and trace- norm 

We make frequent use of two matrix norms [26( in this 
work. The first is the sup-operator norm 



|L4||oo:= sup \\A\v)\\=mj>J{v\AA\v)/y/Wj 

\\\v)\\=l \v) V 

(Al) 

The sup-operator norm of A is the largest eigenvalue of 
\A\ :— VA^A, i.e., the largest singular value of A. Since 
we use it often we denote ||-A||oo f° r simplicity by \\A\\, 
and context should make it clear whether we are taking 
the norm of an operator or simply the Euclidean norm 
lll u )ll := V ( v \ v ) °f a vector \v). Note that if A is nor- 
mal (A^A = AAJ, it can be unitarily diagonalized, so 
that A^A — VD* DV^ where V is unitary and D is the 
diagonal matrix of eigenvalues of A), the largest singu- 
lar value coincides with the largest absolute value of the 
eigenvalues of A, i.e., \\A\\ = supm^^ |(u|A|«)|. 
The trace-norm 

||A||i :=tn/ZU (A2) 

is the sum of the eigenvalues of |j4|, i.e., the sum of the 
singular values of A. Therefore ||A|| < ||-A||x- Both norms 
are unitarily invariant (|| VM.W|| u i = ||-A||ui f° r any pair 
of unitaries V and W) and therefore submultiplicative 
(||AB|| ui < ||A|| ui ||fl|| ui ) In this work we make 

frequent use of both properties. In addition unitarily 
invariant norms are invariant under Hermitian conjuga- 
tion, i.e., ||.A|| u i = ||^4 T || u i. This follows from the singular 
value decomposition: A = VSW^ , where V and W are 
unitaries and £ is the diagonal matrix of singular values 
of A. Since the singular values are all positive we have 
At = WZV^ and hence ||At|| ui = ||S|| ui = ||A|| ui . 



3. Trace-norm distance and fidelity 

The trace-norm distance between two mixed states de- 
scribed by the density operators p\ and p 2 is defined as 



Z?[ Pl ,p 2 ] := -p 2 |U. (A3) 

It is bounded between and 1, vanishes if and only if 
Pi = P2 and is 1 if and only if p\ are p 2 are orthogonal, 
i.e., tr(pip 2 ) = 0. 



The trace-norm distance is a standard and useful mea- 
sure of distinguishability between states. The reason is 
this: Assume that we perform a generalized measure- 
ment (POVM - positive operator valued measurement) 
E with corresponding measurement operators {Ei} satis- 
fying the normalization condition Ei = I. The mea- 
surement outcomes are described by the the measure- 
ment probabilities pi — tr[pi_E!j] and q t — trf^-E^]- The 
Kolmogorov distance between the two probability distri- 
butions produced by these measurements is KE(p,q) = 
3 Tli \Pi ~ an d it can be shown that D[pi 1 p 2 ] = 
max^i^^;, i.e., the trace-norm distance equals the max- 
imum over all possible generalized measurements of the 
Kolmogorov distance between the probability distribu- 
tions resulting from measuring pi and pi [13] • The trace- 
norm distance is related to the Uhlman fidelity 



F[pi,p 2 ] := Wy/plVP^h = toy/y/pip2Vpi (A4) 

via 

1 - D < F < a/1 - D 2 , (A5) 
so that one bounds the other |28| |. 



Appendix B: Proof of the vanishing orders in UDD 

The paper by Yang and Liu [l(| sketches a proof of 
the universality of UDD in the interaction picture. In 
this appendix we provide the first comprehensive proof 
of the Yang-Liu universality result. Our proof is done in 
the toggling frame rather than the interaction picture. 

We shall prove that all powers n < N vanish in the ex- 
pansion of the time evolution operator in (|6|) which have 
an odd number of o~ z B z in Q3, thus also an odd number 
of f z in F3. This is equivalent to showing that the first 
N powers in T of B_(T) vanish, i.e., that dephasing oc- 
curs only in order T N+1 or higher. Henceforth we use 
the shorthand N := N + 1. 

The substitution s = sin 2 (8/2) suggests itself based on 
the UDD choice for the {Sj} in ([1]), because it renders 

f a (0) := f a (sm 2 (e/2)) (Bl) 

particularly simple if the {Sj} are chosen according to 
Eq. Q: 

W) = (B2) 

holds for 6 £ {jir/N, (j + l)n/N) with j G {0, . . . , N}. 
For simplicity, we will omit the tilde on the functions f a 
from now on because only the argument 9 will appear 
henceforth. 

Since f(9) enters the nested integrals only with an 
argument in [0, n] it does not matter what we assume 
about f(9) outside the limits of these integral, and we 
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release the constraint on j, allowing j G Z. The func- 
tion f z (9) then becomes an odd function with antiperiod 
ir/N. Thus its Fourier series 



f z {9) =Y,C2k + ism((2k + l)N9) 



(B3) 



k=0 



contains only harmonics sm(rNO) with r an odd integer. 
The precise coefficients C2k+i do not matter, a fact which 
can be exploited for other purposes, e.g., to deal with 
pulses of finite duration 25 1. 

Under the substitution s = sin 2 (8/2) the infinitesimal 
element ds becomes ds — !• \ sm(9)d6, converting (|6b[) into 



Fx = 



(B4) 



sm(0 n )d0 n f a7l (0 n ) / sisx(0n-i)d0n-ifa n ^ 1 (O n -i) 
o Jo 



sin(0 2 )d0 2 / Q2 (0 2 ) / sm{e 1 )de 1 f ai (e 1 ), 

o Jo 



where we absorbed the h factors coming from the in- 
finitesimal elements into the coefficients c 2 fe+i. 

What happens if we perform the successive integra- 
tions in Eq. (|B4|) ? Replacing f z (9) by its Fourier se- 
ries (|B3[) we deal with integrands which are products of 
trigonometric functions. The substitution gave rise to 
the factor sin 9. The Fourier series gives rise to additional 
factors sm(r o N0), where r D is an odd integer. Recall the 
elementary trigonometric identities 

sin a sin 6 = — [cos (a — b) — cos (a + b)] , (B5a) 
cos a sin b — — [sin (a + b) — sin (a — b)] , (B5b) 
cos a cos b = — [cos (a + b) + cos (a — b)] . (B5c) 

Using these, the most general trigonometric factor to oc- 
cur in the course of the integrations in Eq. (|B4|) can be 
written as either sm[(q + rN)0] or cos[(q + rN)9] where r 
and q are integers. Since we are only concerned with val- 
ues of n such that n < N, the absolute value of q always 
remains below N, so that the representation of the inte- 
ger factor (q + rN) in the arguments of the trigonometric 
functions is unique. 

We now consider a complete set of four different cases 
which can occur in the course of the evaluation of each 
F$. The first two cases are associated with the occur- 
rence of fo — 1 in one or more of the nested integrals. 
Suppose for concreteness that this happens in nested in- 
tegral number j. Then the factor sm(r N9j) does not 
occur in this integral, since this factor arises exclusively 
due to the presence of f z {9j)- The two cases are now 
distinguished by whether a summand in the integrand of 
this jth integral, after a complete expansion of trigono- 
metric products excluding the sin(9j) term, into sums us- 
ing Eqs. (IB5a[) - (|B5cl) . involves the factor cos[(g + rN)9] 



(whence we call the integrand "cosine- type" ) or the fac- 
tor sin[(g + rN)9] (whence we call the integrand "sine- 
type"). The third and fourth cases are associated with 
the occurrence of f z in integrand number j. Then the 
factor sin(r N9j) does occur in this integrand, and again 
we distinguish two cases according to the presence of 
cos[(g + rN)9] ( "cosine- type" ) or sin[(q + rN)9] ("sine- 
type") arising from a complete expansion of trigono- 
metric products excluding the sin(9j) term and also the 
s'm(r N9j) term. Here, then, are the four cases in detail: 



(i) Assume that one of the nested integrals contains 
fo = 1 and the factor cos[(q + rN)9]. As we shall 
see in item 4) below this case occurs when r is even. 
Then this integral reads 



2 / cos[(<7 + rN)9] sm{0)d9 = (B6) 
cos[(-l + q + rN)9] cos[(l + q + rN)0] 



-1 + q + rN 



I + q+rN 



In writing Eq. (|B6[) we have assumed that the de- 
nominators do not vanish. The denominators may 
in fact vanish because r may be zero. When r = 
the case \q\ = 1 is special and yields 

2 J cos(±0)sin(9)d9 = -i C os(20). (B7) 



The important point is that both Eq. (IB6|) and (|B7|) 
have only cosine terms on the right hand side. 



(ii) Assume that one of the nested integrals contains 
fo = 1 and the factor sin[(q + rN)0]. As we shall 
see in item 4) below this case occurs when r is odd. 
Then this integral reads 



2 / sin[(<? + rN)9] sva(9)d0 = (B8) 
sin[(-l + q + rN)0] sin[(l + q + rN)0] 



-1 



rN 



1+q+rN 



No denominator can vanish because as we shall see 
in item 2) below, \q\ < N. The important point 
here is that Eq. (|B8[) has only sine terms on the 
right hand side. 



(iii) Assume that one of the nested integrals contains 
f z and the factor cos[(g + rN)9]. As we shall see 
in item 4) below this case occurs when r is even. 
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Then this integral reads 



4 / cos[(<? + rN)9] sm(r o N0) sm(9)d9 

sin[(-l + q + {r + r Q )N)9] 
-1 + q + (r + r a )N 
sin[(l + q + (r + r a )N)9] 

1 + q + [r + r a )N 
sin[(-l + g + (r - r o )N)0] 

-l + q+(r- r a )N 
sin[(l + g + (r - r )N)9] 
l + q+(r- r Q )N 



(B9) 



Since r ± r c is odd none of the denominators can 
vanish as long as \q\ < N. Again, the important 
point here is that Eq. (|B9|) has only sine terms on 
the right hand side. 

(iv) Finally, assume that one of the nested integrals con- 
tains f z and the factor sin[(g + rN)9\. As we shall 
see in item 4) below this case occurs when r is odd. 
Then this integral reads 



sin[(g + rN)9] sin(r N9) sin(0)d0 

cos[(-l + q + (r ~ r o )N)0] 
-l + q+(r- r a )N 
cos[(l + q + (r - r ) N)9] 

1 + q + (r - r a )N 
cos[(-l + q+(r + r )N)9] 

-1 + q + (r + r )N 
cos[(l + q + (r + r )N)9] 
l + q + (r + r )N 



(BIO) 



In writing Eq. (|B10I) we have assumed that the 
denominators do not vanish. A denominator can 
vanish only when |g| = 1, which leads to the two 
special cases r = ±r G . In analyzing these two cases 
we can assume without loss of generality that q = 1 
and r = r . Otherwise we multiply the argument 
of the first and/or the second sine- function in the 
integrand by — 1. This yields 

-I I sin[(l + rN)9] sm(rN6) sin(9)d9 = (Bll) 
cos(26>) cos(2rN9) cos[(2 + 2rN)0] 



2rN 



2rN 



The important point here is that in Eq. (|B10|) only 
cosine terms appear on the right hand side. 

The number of possible terms proliferates in the course 
of the successive integrations. Therefore, in the sequel we 
discuss only the common features of the resulting sum- 
mands. It is always understood that sums with varying 
sets of q and r are considered. We present a series of 



observations which leads to the desired proof of the can- 
cellation of the first N powers in T of B_(T). The key 
to doing so will be to show that after integrating with 
an odd number of f z factors we always end up with a 
sine-type integrand. 

1) Recall that we call an integrand summand "cosine- 
type" or "sine-type" if after complete expansion 
of all trigonometric products, excluding the sin(0) 
term arising from the change of variables and of 
the f z term if it is there, the trigonometric factor is 
cos[(q + rN)9] or sm[(q+rN)9], respectively. Cases 
(i) and (iii) above are cosine-type, while cases (ii) 
and (iv) are sine-type. It is clear that the first inte- 
grand in (IB4[) is cosine- type with r = and q = 0. 

2) We track which values of q may occur in each inte- 
gration. The first integrand in (|B4[) is either sin(0i) 
or sin(0i) sin(rA^i), i.e., it has q — 0, so that this 
is our starting point. It follows from Eqs. (IB6|) - 
(jBllI) that each integration increments the possi- 
ble maximum of \q\ by unity. The highest power in 
T studied is T N so that there are n < N integra- 
tions. This implies that the final value qa ns d before 
the very last integration obeys | Qfinai | < N. 

3) We track whether even or odd values of r occur at 
each integration. As mentioned in item 1), r = 
holds in the first integration, so that our starting 
point is an even value. Each integration involving 
fo leaves r unchanged. Each integration involving 
f z adds or substracts r a , so that r changes from 
even to odd or vice versa. If we combine this with 
the results of cases (i)-(iv) [Eqs. (|E5|)-(|BTT|) this re- 
veals the input-output table in Eq. (|B12[) . The in- 
tegrands (i.e., the inputs) are indicated by the case 
number in the table entries, while the values of the 
integrals (i.e., the outputs) are the types indicated 
in corresponding table entries. Also indicated is 
the transformation undergone by r from input to 
output. 



X 


cosine- type 


sine-type 


fo 


case (i): 


case (ii): 




cosine-type, r M> r 


sine- type, r > r 


fz 


case (iii): 


case (iv): 




sine-type, r t— > r ± r 


cosine-type, r t— > r ± r D 



(B12) 

4) Consider the output of the table as the input into 
the next integration and focus on the f z row. Note 
that cases (iii) and (iv) alternate along with a 
change in parity of r, i.e., cosine- type changes into 
sine-type and vice versa, while odd r changes to 
even r and vice versa. Therefore if we start with a 
cosine-type integrand and perform an odd number 
°f fz integrations, we will end up with a sine-type 
output and a change in parity of r. For the same 
reason, since the first integrand in (|B4I) is cosine- 
type with even r, and case (i) can only be arrived at 
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after an even number of f z integrations (the num- 
ber of /o integrations is arbitrary), case (i) always 
involves even r. Repeating this reasoning explains 
why case (iii) also has even r, while cases (ii) and 
(iv) have odd r. 

5) Dephasing results only from the terms which com- 
prise an odd number of f z integrations. Consider- 
ing that as noted in item 1) we start from a cosine- 
type integral and with r = 0, it follows from item 
4) that the last integration provides a sine-type re- 
sult. This integral can therefore be written as a 
sum over terms all of which are of the form 



sin[((7fi na i + 1 + r fina iiV)0] 



0. 



(B13) 



We used the definition of b a p{T) [Eq. ([26]) ] in the last 
equality. Next, we use the triangle inequality, and finally 
the unitary invariance of the trace norm along with the 
normalization of \ib) 



2D[p s (T),pUT)} 

<|6 ++ (T)-1||||^|||i + |6+ 
+ |&_ + (T)|||a z |V>)<V|||i + |&- 
= |6 ++ (T)-1| + |6 + _(T)| + |6. 



_(T)| 

. + (T)| + |6__(T)|. 



(C2) 



Recall that the operator i?_ in Eq. (fTTj) comprises all 
the terms with odd number of o~ z ® B z . Hence we have 
proven that the first N powers in T of B_(T) vanish, i.e., 
B_{T) = 0(T N+1 ). This is what we set out to show and 
concludes the derivation. 

A remark concerning the result of Yang and Liu ob- 
tained in the interaction picture (lOj is in order. They 
showed that exp(iT Bq)U (T) comprises only odd pow- 
ers in a z which are of order T N+1 or higher. Since 
exp(±iTB ) does not contain any term proportional to 
a z the Yang-Liu result implies our result and vice versa. 

For time-dependent Hamiltonian, the proof in the in- 
teraction picture Q is more convenient because powers 
in time occur anyway. Therefore we stress that the state- 
ment that only odd powers in u z of order T N+1 or higher 
occur is independent of the choice of reference frame, i.e., 
the description in the interaction picture or in the tog- 
gling frame. 



Appendix D: Correlation function inequality 



We prove the correlation function inequality (|32|) . 
Consider the spectral decomposition of the bath density 
operator: pb = where {\i}} are normalized 

eigenstates, A-j > are the eigenvalues, and ^ A-j = 1. 
Defining := Q\i) and \v!j) := (Q')^\i), we have in this 
eigenbasis of pb- 



Appendix C: Distance bound calculation 

We prove the trace- norm distance bound Eq. ([27)) : 



2D[p s (T),p° s (T)] 
= ||tr B [psb(T)-p° S b(T)} ||i 
= ||tr B [U(T)p° SB (0)U\T)] - tr B [pg B (0)] || x 
= litre [(/ ® B+{T) + a z ® B _{T)){\i>) ® p B 

(/®B|(T) + o-,®Bl(T))] -WWWi 

= \\(b ++ (T) - l)\i;)^\+b + .(TM{^\a z 
+ b-+(T)a z \^\+b—(T)a z \^\a g \\ 1 . (CI) 



|tr [Q P bQ'\ I = |tr [Q'Qpb] | = | ^{i\Q'Q\i)\ 

i 

i i i 

<Eiiq'iii^ii a ^iiq'iiiiqii ( D1 ) 



where we used the triangle inequality, followed by 
the Cauchy-Schwartz inequality, and then the bounds 
Ilk>|| = ||Q|*)|| < IIQH and || - \\{Q'm\\ < 
ll(Q')^l = IIQ'IIj which follow from the definition and 
properties of the sup-operator norm (see Appendix |A"]) . 
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